Critical fields and devil's staircase in superconducting ladders 
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We have determined the ground state for both a ladder 
array of Josephson junctions and a ladder of thin supercon- 
ducting wires. We find that the repulsive interaction between 
vortices falls off exponentially with separation. The fact that 
the interaction is short-range leads to novel phenomena. The 
ground state vortex density exhibits a complete devil's stair- 
case as the applied magnetic field is increased, each step pro- 
ducing a pair of metal-insulator transitions. The critical fields 
in the staircase are all calculated analytically and depend only 
on the connectivity of the ladder and the area of the elemen- 
tary plaquette. In particular the normal square ladder con- 
tains no vortices at all until the flux per plaquette reaches 

2V3' 



It is now relatively easy to fabricate complex structures 
containing many Josephson junctions. Two-dimensional 
arrays have long been of interest. In a transverse mag- 
netic field, both square ||] and triangular [|[j4) 2D arrays 
show very rich structure in resistance, impedance and in- 
verse sheet inductance measurements when the flux per 
plaquette is close to a simple rational fraction (at which 
commensurate states form). Such arrays appear to be 
well described by a frustrated XY model |5Hf!| although 
special attention has only been given |jMto the sim- 
plest rational flux values: / = |, | and g7 Vortex dy- 
namics in 2D arrays display many anomalous features 



[15,Q. Interesting self- inductance |L7|,[18| and charging 
[ 19 1 effects arc also observed. 

Analogous phenomena are observed in superconduct- 
ing weakly coupled wire networks made in thin amor- 
phous niobium silicon films pOLpll . Furthermore Monte- 



Carlo calculations |22| on a Coulomb gas model of a 2D 
network have indicated the existence of two melting tran- 
sitions in the vortex sub-system. 

Recently there has been much interest in ladder arrays 
of Josephson junctions. Experimental work p^j on wide 
ladders has found a rich structure in the dependence of 
the longitudinal resistance on transverse magnetic field 
and this has been interpreted in terms of metal-insulator 
transitions. In numerical calculations ]24| it has been no- 
ticed that the normal square Josephson ladder exhibits 
a critical field below which the ground state contains no 
vortices while Monte-Carlo simulations have shown, for 
the case / = k, that the ladder belongs |25| to the uni- 



versality class of the Ising model. 

A quasi-classical analytical method for investigating 



Josephson networks has been successfully tested pq] 
against full Monte-Carlo simulations. In this paper we 
use the same method to completely determine the ground 
state, as a function of applied magnetic field, of both the 
Josephson ladder shown in Fig. la and of the analogous 
superconducting wire network shown in Fig. lb. We find 
that the vortex-vortex interaction falls off exponentially 
with separation, in contrast to the long-range logarith- 
mic interaction between vortices in 2D and 3D arrays, 
and that at low temperatures the ladder is described by 
an Ising model. The short-range nature of the interaction 
gives rise to a geometry dependent critical field H g which 
depends solely on the geometrical structure of the ladder 
and plays the same role as H c \ in bulk superconductors. 
Above H g the vortex density exhibits a complete devil's 
staircase as the applied field is increased. To our knowl- 
edge, this is the first example of a real physical system 
exhibiting a non-trivial complete devil's staircase. 

First we study the Josephson ladder and then show 
that the analogous ladder of superconducting wires be- 
haves in the same way. A network of Josephson junctions 
can be described by the XY- model [E7| : 
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where the sum is over all junctions, I c is the critical 
current and $0 = is the flux quantum. The gauge- 
invariant phase difference 4n across the ith junction is 
defined by 
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where A8i is the phase difference and the line integral of 
the vector potential A is taken across the junction. The 
current through the ith junction is given by 



Ii = I c sin < 



(3) 



Consider a ladder (Fig. la), which we call the Joseph- 
son ladder, consisting of a linear series of identical pla- 
quettes. Let there be 2u unshared junctions in each pla- 
quette and c junctions in each region which is in common 
between two plaquettes. Let the applied magnetic flux 
per plaquette be f&o- The sum of the gauge invariant 
phase differences around any plaquette is then quantised 
according to: 
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where rij is the vortex number and 2$ is the current 
which, if circulating around a plaquette, would give an 
induced flux of one flux quantum $ - The last term gives 
a good approximation to the induced flux if, for example, 
the junctions are uniformly spaced around a rectangle. 

If u and if the number of junctions in each plaque- 
tte is large then the phase <f> dropped across each junc- 
tion will be small enough that we have sin <fi ss <fi and so 
Eqs. (§,§ give us 



£ hR = n. } - f 

plaquette j 



where the 'resistance' R is 

r= 1 
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In this approximation the energy of junction i becomes 
/ l 2 $o/(47r-fc) and so the total energy U is 



U = K 



all junctions 



(7) 



where n = <5>o/(4ttI c ). Since the currents are linear in / 
the energy U must be quadratic in /. 

Now we consider a similar ladder of superconduct- 
ing wires as shown in Fig. lb. ua and ca are the 
lengths of the wires, u and c being integers. When 
the true superconducting state is established the su- 
perconducting order parameter takes the form ijj = 

V^exp ji (d - |£ J A.dl^j. If the width of each wire 
is small compared with the London penetration depth 
then vortices cannot penetrate the wires; if in addition 
the width is small compared with the length then the 
current density J will be constant: 



J 
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where S is the wire cross-sectional area and / is the cur- 
rent carried by the wire. Integrating around a plaquette 
we again obtain Eq. @ but now we have 



R = 



1 
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The total energy of the ladder is likewise given by Eq. (Q) 
except that now k — ml (4e 2 |Vo| 2 S). 

Thus Eqs. (||) and (Q) determine the solution of the 
currents and total energy in both the Josephson ladder 
and a ladder of long thin superconducting wires (and also, 
of course, the solution of the currents and power dissipa- 
tion in a similar network of resistors and batteries). This 
method has already been tested J2q] against Monte-Carlo 
calculations for a number of small Josephson networks. 

Solving Eqs. (^) and ([?]) for a ladder of infinite length 
we find that the energy U is 



(10) 
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where ej = rij — /, Uq = k/ (4i? 2 V u 2 + 2uc) and 



, 2u u 
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a depends solely on u/c as shown in Fig. 2; all other 
factors contribute only to the energy scale Uq. Note that 
vortices repel each other (a < 1) and since the energy de- 
creases exponentially with separation the vortex-vortex 
interaction is a short-range interaction of range — 1/ In a 
plaquettes. 

We now focus on finding the lowest energy state for 
any given value of /. For / in the range < / < 1 we 
need only consider two values of rij , namely rij — and 
rij = 1. Hubbard p8| has shown that for particles (in 
our case vortices, i.e. plaquettes with rij = 1) interact- 
ing with a convex repulsive interaction the lowest energy 
state of p particles in a one-dimensional chain on q sites 
is a generalised Wigner lattice consisting of a sequence of 
two spacings, the spacings being the two integers which 
bracket q/p (for example, when q/p = 7/3 the two spac- 
ings will be 2 and 3). Fig. 3 shows the energies of these 
generalised Wigner lattices for various vortex densities 
p = p/q as calculated by direct solution of Eq. (||) . Imme- 
diately we see that there is a critical flux f c below which 
the empty ladder (i.e. p = 0) is the lowest energy state. 
A more accurate numerical solution was performed for 
the same Josephson ladder, fa in Eq. ([|) being replaced 
by sin -1 (/;// c ) instead of Ii/I c ; this more accurate calcu- 
lation gave the transitions from one generalised Wigner 
lattice to another as occurring at much the same values 
of /, differing by no more than 0.01 from the values seen 
in Fig. 3. A ladder of larger u and c would have shown 
better agreement. 

We can calculate the critical flux f c very simply. The 
energy increase 5U when one vortex is added to an empty 
ladder is 



su _ 1 2 i + t* 

2U 1 - a 



f 



(9) Thus 8U > for / < f c where 
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i.e. for / < f c the lowest energy state is the empty 
ladder. Furthermore we see that, as shown in Fig 2, the 
critical flux f c depends solely on the ladder ratio u/c; 
it is independent of all other parameters. Eq. ( p^| ) is 
consistent with the result obtained previously 0] for the 
square ladder with nearest neighbour interactions only. 

It may look as though Fig 3 not only tells us the ground 
state but also the low energy excited states. However 
this is not the case. Excited states are not necessarily 
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generalised Wigner lattices; in fact since the interaction 
is short range the lowest energy excitation will normally 
be a domain wall. 

Since, for low energy states, there are only two possible 
values of e (i.e. s = — / and e = 1 — /) Eq. ( |i~0| ) is 
reminiscent of the Ising model Hamiltonian. Let 
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Re-writing Eq. ( J10| ) in terms of these Ising variables we 
obtain for a ladder of N plaquettes (valid in the limit 

N -> oo) 
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This is an anti-ferromagnetic Ising model Hamiltonian 
in which \ — f plays the role of an applied field, the 
interaction between spins falling off exponentially with 
separation. The fact that at low temperatures the ladder 
is described by an Ising model Hamiltonian explains the 
Monte-Carlo result Jl3| (determined only for the special 
case / = 1/2) that the Josephson ladder belongs to the 
universality class of the Ising Model. 

We now explore the implications of Eq. (|l5|). Pub- 
lished literature p9|-]3l| gives no complete solution of the 
thermodynamics of this kind of Ising model at both non- 
zero field and non-zero temperature. It has been shown 
p2| however that the ground state of a one-dimensional 
Ising model with an anti-ferromagnetic interaction ex- 
hibits a complete devil's staircase [[53| if the interaction 
is convex (i.e. decreases with separation with positive 
2nd derivative); the average spin (s) per site is always a 
rational fraction, all rational fractions between and -1 
being visited as the field (in our case \ ~ f) is increased 
from zero to the critical field ( in our case ^ — f c ) above 
which the spins are all aligned opposite to the field (in 
our case the ladder contains no vortices). 

Krantz et al |54| have calculated the precise nature 
of the expected devil's staircase for the case of an anti- 
ferromagnetic Ising model with an exponential interac- 
tion between spins. Using Eq. (15) their results may now 
be applied to the superconducting ladder. Fig. 3 shows 
the predicted devil's staircase for the 'square' ladder (i.e. 
u = c = 1); note that the locations of the transition 
points are seen to agree with the crossing points of the 
energy curves. A ladder with u < c would give a stair- 
case in which steps corresponding to larger values of the 
denominator q would be more easily resolved. 

In conclusion, we have shown that in the Josephson 
ladder and in the superconducting wire ladder the in- 
teraction between vortices falls exponentially with sepa- 
ration; this is in contrast to the 2D array, where the 



logarithmic interaction between vortices leads to the 
Kosterlitz-Thouless phase transition j3j|. At low tem- 
peratures the ladder is described by a ID Ising model 
(Eq. (15)) and so no phase transition is expected |3(| ex- 
cept at zero temperature. This also leads to the ground 
state vortex density exhibiting a complete devil's staircase 
as the applied magnetic field is increased] at each step in 
the staircase one extra vortex penetrates the ladder. To 
our knowledge this is the first example of a physical sys- 
tem exhibiting a non-trivial devil's staircase. The fact 
that the devil's staircase occurs in both the Josephson 
ladder and the wire ladder suggests that it is a very gen- 
eral property of superconducting ladder structures. 

In 2D networks and in wide multi-leg ladders [p3[ com- 
mensurate states have been observed centred on ratio- 
nal flux values; note that our results suggest that this is 
not the case for the ID ladder, e.g. a vortex density of 
p = 1/5 doesn't occur at / = 1/5 but in a small region 
close to / = 0.3. Likewise, the critical field f c has no ap- 
parent analogue in 2D. These important differences arise 
because in the simple ID ladder the interaction is short- 
range (exponential) whereas in 2D it is logarithmic. It is 
interesting that in the long range limit (i.e. small u/c) 
the properties of the ID ladder are reminiscent of the 2D 
ladder: i.e. f c — > and the ground state vortex density 
p = p/q^ f. 

We call for experiments to detect the predicted devil's 
staircase phenomenon which should be observable in low 
temperature mutual inductance and resistivity measure- 
ments. Each step in the staircase should give rise to a 
pair of metal-insulator transitions p3| : close to a step 
the disorder will cause high resistance while far away 
from any step the ordered array of vortices will cause 
low resistance. 

From a practical point of view perhaps the most im- 
portant result is that both kinds of ladder exhibit a 
geometry- dependent critical field H g , analagous to H&, 
below which all vortices are expelled. Applying Eqs. (|Tl| ) 
and (|l3| ) one can design a ladder of any desired H g = 
fc<&o/A by appropriately choosing the ladder ratio u/c 
and the area A of the elementary plaquette. 

We would like to thank John Samson and Daniil 
Khomskii for illuminating discussions. 
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FIG. 3. The solid curves show the dependence of en- 
ergy per plaquette U/(NUo) on flux per plaquette / for 
the 'square' ladder, i.e. u = c = 1; each curve repre- 
sents a generalised Wigner lattices of a particular vortex den- 
sity p — p/q. The energies are calculated for a 512 plaquette 
ladder by direct solution of Eqs. (|E]) for the case 7$ = oo. 
The broken line shows the predicted devil's staircase in the 
ground state vortex density; each visible step contains an in- 
finite number of unresolved steps. Note the critical flux at 
/ = f c = l/(2v^) « 0.289. 



FIG. 1. Various ladder structures: (a) the general Joseph- 
son ladder, (b) the wire ladder. 



FIG. 2. The dependence of f c and a on the ladder ratio 
u/c. 
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